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ABSTRACT 

We  construct  compactly  supported  wavelet  bases  satisfying  homogeneous 
boundary  conditions  on  the  interval  [0,1].  The  maximum  features  of  multires¬ 
olution  analysis  on  the  line  are  retained,  including  polynomial  approximation 
and  tree  algorithms.  The  case  of  i^oliO,  1])  is  detailed,  and  numerical  values, 
required  for  the  implementation,  are  provided  for  the  Neumann  and  Dirichlet 
boundary  conditions. 
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INTRODUCTION 

Wavelet  bases  are  often  presented  as  a  powerful  tool  to  perform  the  ap¬ 
proximation  and  the  numerical  resolution  of  partial  differential  equations. 
Indeed,  thanks  to  zero  moment  and  localization  properties,  wavelet  spaces 
are  self-adapted  to  the  solution  and  therefore  may  allow  fast  and  accurate 
resolution.  In  the  last  few  years  different  algorithms  have  been  successfully 
tested  on  linear  and  non-linear  equations  [9]  [12].  Nevertheless,  most  of  the 
addressed  problems  where  posed  in  the  periodic  framework  which  circum¬ 
vents  the  difficulties  generated  by  general  boundary  conditions,  but  which 
makes  nearly  impossible  the  treatments  of  real  problems. 

It  is  known  [3]  that  under  very  general  hypotheses,  one  can  only  consider 
homogeneous  conditions.  Therefore,  we  are  driven  to  the  construction  of  ap¬ 
proximation  spaces  for  homogeneous  functional  spaces  on  the  interval.  For 
example,  for  problem  involving  homogeneous  Dirichlet  conditions  on  [0, 1], 
i.e  tt(0)  =  w(l)  =  0,  the  solution  u  could  be  reached  in  the  Sobolev  space 
ffo([0, 1])  =  {u  6  l]),u(0)  =  u(l)  =  0}  and  this  leads  to  the  construc¬ 

tion  of  approximation  spaces  for  Hq{[0,  1]).  Various  constructions  have  been 
proposed  (see  Auscher  [2],  in  [4]  and  references  in  it)  but  are  not,  according  to 
their  authors,  numerically  tractable.  The  proposed  construction  follows  the 
one  of  Cohen  et  al.  ([1]).  After  a  short  recall  of  this  preliminary  construction 
we  introduce  and  analyze  the  homogeneous  space  construction.  Numerical 
details  are  provided. 


I  THE  PRELIMINARY  CONSTRUCTION 
OF  COHEN,  DAUBECHIES  AND  VIAL 


This  construction  of  wavelets  on  the  interval,  ([!]),  is  derived  from  the  com¬ 
pactly  supported  wavelet  multi  resolution  analysis  on  the  line  introduced  by 
I.  Daubechies  [5]. 

In  the  case  of  L^{1R)  this  multi  resolution  analysis  is  classically  given  by 
a  sequence  of  closed  subspaces  Vj  satisfying: 

i)  ...  C  V-i  C  Vo  C  Vi  C  V2  C  ...L\R),  CijezVj  =  {0}  and  = 

L^JR) 
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Moreover,  each  Vj  is  spanned  by  the  translations  of  the  dilated  version  of  a 
fixed  function,  the  scaling  function  <f>,  i.e  Vj  =  span{2^^'^(j>{2K  —  fe),  k  £  Z}. 
Here,  the  family  {2^^'^4>{2^ .  -  k),  k  e  Z}  is  orthonormal  and  (f>  is  on  the 
one  hand  compactly  supported  and  on  the  other  hand  such  that  its  Fourier 
transform  ^  satisfies  the  Strang  and  Fix  approximation  rules  of  order  N-1 

<j>^'^\2kir)  =  0,  A:€^\{0},  n  =  0, ...,  iV  -  1.  (1) 

One  consequence  of  (1)  is  that  the  family,  {<!){.-  k),  k  €  Z},  can 
reproduce  locally  the  polynomials  of  degree  at  most  —  1. 

The  support  of  ^  is  the  interval  [-A^  +  1,  A^]  and  the  regularity  of  (f>  is 
asymptotically  [5].  Moreover,  (j)  is  solution  of  the  following  scaling 

equation: 


hk(l){2x  -  k).  (2) 

jk=-iV+l 

The  detail  spaces  Wj  are  defined  as  the  orthogonal  complements  of  Vj  in 
^•+1,  i.e. 


Wi  =  n 


(3) 


and,  thanks  to  i) 

(JV  =  L\K). 

j^Z 

The  essential  feature  of  multi  resolution  analysis  (see  Y.  Meyer  [11])  is  that 
3V’  such  that  Vj  €  2Z 

Wj  =  span{2^^^‘ij}{2K  —  k),  k  £  Z}. 

Again,  the  family  {2^l‘^ij){2K  —  k),  k  £  2Z}  is  orthonormal.  The  function 
V’  is  here  a  compactly  supported  wavelet  and  is  obtained  from  the  following 
detail  equation 


^{x)=  Yu  Sk(i>{2x-k)  (4) 

fc=-JV+l 
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Moreover  $upp{il>)  =  supp{<f>)  and  V’  has  the  same  regularity  as  (j).  In  addition, 
because  of  the  approximation  properties  of  Vq  and  the  definition  of  Wo, 
has  got  N  vanishing  moments,  i.e: 

jx^‘ip{x)  =  0  I  =  0,...,N  —  1. 

Finally,  the  family  (2-' a:  —  A:)}  is  an  unconditional  basis  for 

various  functional  spaces  such  as  Holder  spaces  C^{1R.)  or  Sobolev  spaces 

In  [1],  the  goal  of  I.  Daubechies  et  al.  was  to  construct  a  family  of  wavelet 
basis  on  the  interval  [0,1]  able  to  characterize  L'^{[0, 1]),  1])  or  (^*([0, 1]) 

while  preserving  the  most  attractive  properties  of  multi  resolution  analysis 
of  L‘^{]R),  despite  the  lack  of  shift  invariance  of  L^{[0, 1])  (this  is  not  the  case 
for  the  constructions  of  P.  Auscher  ( in  [4])). 

We  give  in  the  following  paragraph  an  outline  of  the  construction  but  the 
reader  should  refer  to  [1]  for  details. 

The  construction  is  performed  in  two  steps  as  follows. 

The  first  step  consists  in  defining  suitable  subspaces  of  T^([0,  Ij)  from 
a  basis  essentially  constructed  from  the  translated  versions  of  a  rescaled 
function,  while  the  second  step  consists  in  the  construction  of  the  detail 
spaces  with  the  same  requirement. 

More  precisely,  in  the  first  step,  Vj([0, 1])  is  constructed  as  follows: 

Thanks  to  the  compact  support  of  for  large  enough  values  of  j  and 
k  =  N,  ...,2^  —  —  1,  the  support  of  the  functions  (l>{2^x  —  k)  is  included  in 

[0,1].  Therefore,  the  corresponding  functions  may  be  used  as  the  interior  basis 
functions  of  I^([0, 1])  and  the  set  =  {(^(2-^x  —  A;),  k  =  N,  —  N  —1}  (/ 
stands  for  interior)  is  then  defined.  To  fully  define  Vj([0, 1]),  N  edge  functions 
are  added  at  each  boundary  of  [0,1]  to  complete  the  basis  $/.  These  two 
families  of  N  functions,  ^e,o  =  k  =  0,  ...,7V— 1}  and  $£,1  =  k  = 

2^  —  TV,  ...,2^  —  1}  are  constructed  to  have  minimal  support  and,  such  that 
the  order  of  approximation,  (TV),  related  to  the  interior  functions  is  kept. 
In  other  words,  all  polynomials  of  degree  less  than  TV  —  1  should  be  locally 
expandable  as  a  linear  combination  of  the  basis  functions  of  Vj([0, 1]).  Let  us 

^We  remind  that  ,  s  €  jR,  /  belongs  to  H^{T)  if  and  only  if  l/(”)P(l  +  ”^)*  < 

+00  and  that  for  0  <  a  <  1  ,  /  €  C"  if  and  only  if  \f(x  +  h)  —  f(x)\  <  C\h\°‘  for  every 
X,  h  in  JR,  the  constant  C  not  depending  on  x  and  h. 
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recall  the  construction  of  ^e,o  (the  same  applies  to  The  edge  functions 

^E,o  are  defined  for  k  —  Q,...,N  —  1  as  the  restriction  to  [0,1]  of  a  specific 
linear  combination  of  the  family  {<f>{x  —  ^)such  that  0  €  supp[<f>{x  —  fc)}. 
More  precisely, 

2N-2 

=  2-'^^  Y.  (J  <^(2^2:  +  n  -  iV  +  1)  X[o,i]  (5) 


The  supports  supp{<p]^j^)  are  staggered,  i.e,  supp{ip°k)  =  [0. 


27V-1-A:i 


and  (fill, 


is  a  polynomial  of  degree  k  on  the  interval  [0,^]  (see  Figure  1  for  example). 
Indeed,  —  n)  is  a  polynomial  of  order  k  [1]  and  is  a  polynomial 

in  n  of  degree  k. 

By  construction,  ^e,o  -L  but  $£;,o  is  not  an  orthonormal  family.  An 
orthonormalization  procedure  using  the  Gram-Schmidt  algorithm  is  then  per¬ 
formed.  Starting  from  jv_i  down  to  Q  one  obtains  N  orthonormal  edge 
functions  ^  =  0, ..., iV—  1}  with  staggered  support  [0, iV-h  A:]  and  still, 

2-j]  is  a  polynomial  of  degree  k. 

Finally,  ij([0, 1])  is  then  by  definition  generated  by  the  orthonormal  family 
^jE,o  U  U  $£,1: 

f  A;  =  0,...,iV-l}  ] 

U 

1^([0, 1])  =  span  <  {(t)j,k,  k  =  N,...,2^  -  N  -  1}  >  (6) 

u 

I  fc  =  2^  -7V,...,2^  -l}  J 

with  4>jk  =  2^l‘^4>{2^x  —  k).  One  gets 

^o([0,l])  c  V,o+i([0,l])  C  ...  C  v;([0,l])  C  ...T^([0,1]) 

where  jo  is  chosen  so  that  support{^E,o)  0  support{^E,i)  =  0-  As  <l>jk,  the 
edge  functions  satisfy  a  modified  scaling  equation  (2),  one  writes 


N-l  N+2k 

Y  +  Y  ^l,n4>j+hr> 

n=0  n=N 


k  =  0,...,N-l  (7) 


The  numerical  values  of  the  coefficients  n  =  0,  ...,7V  -f-  2A;;  k  = 

0,...,Ar-l}  and  {hl^,  n  =  3*2^ -N-2k-3,’...,2^ -1-,  k  =  27-iV,  ...,27-1} 


4 


for  the  right  edge,  are  computed  in  [1]. 


The  second  step  of  the  construction  is  the  definition  of  a  suitable  basis 
for  the  usual  wavelet  space  lTj([0, 1])  =  Frti([0,l])n(VS([0.11))"’-  Thanks 
to  (3)  and  to  the  compact  support  of  for  large  values  of  j  ,  the  family 
=  {^(2'' a:  —  fc),  k  ■=  —  N  —  1}  belong  to  Wj([0, 1]).  Wavelets 

of  this  family  constitute  a  first  part  of  the  basis  of  VFj([0, 1])  and  are  called 
the  interior  wavelets.  Since  dim{Wj{[0,l])  =  2^,  N  other  wavelets  at  each 
edge  should  be  added  to  $/.  Again,  we  only  recall  the  construction  at  the 
edge  a;  =  0.  The  complementary  wavelets  are  deduced  from  the  definition  of 
W,([0, 1])  as 

=  vln  k  =  (8) 

n=0 

where  (., .)  stands  for  the  scalar  product  of  T^([0, 1]).  By  construction  they 
are  orthogonal  to  V^  ([0, 1])  and  to  the  interior  wavelets.  Their  supports  are  no 
longer  staggered,  but  an  iterative  process  described  in  [1]  reduces  the  support 
of  to  [0,^^^]  instead  of  The  last  step  of  this  construction 

consists  again  of  a  Gram-Schmidt  orthonormalization.  Starting  from  k  =  0 
up  to  A— 1  one  gets  N  orthonormalized  wavelets  for  the  left  edge,  {tpjk  ,  k  = 
0,  ...,N  —  1}.  These  wavelets  are  known  through  the  coefficients  n  = 

0, ...,  N  +  2k-,  k  =  0, ...,  iV  —  1}  that  occur  in  the  modified  details  equation: 

N-l  N+2k 

=  X]  ^  =  0,  ...,  N  (9) 

^=0  7Z=iV 

W,([0, 1])  is  therefore  entirely  characterized  by 

r  k  =  D,...,N-i}  1 

u 

W,([0,1])  =  span{  k  =  N,...,2^  -N  i  (10) 

U 

I  k  =  2i-N,...,2^-\)  J 

Since 

imii)=v'»ao.ii)©Mo([o,i]) 

3>k 
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one  gets  an  orthonormal  basis  of  L^{[0, 1])  as: 


f  fc  =  o,...,iv-i}  1 

&  =  o,...,i\r-i}  1 

u 

U 

u 

i>jo 

[J 

V 

II 

to 

o 

1 

to 

o 

1 

1 

1 

II 

Remarks: 

As  we  have  said  before,  these  wavelets  bases  are  very  attractive  because 
they  preserve  the  main  features  of  the  whole  line  construction.  More  pre¬ 
cisely,  since  the  edge  functions  are  finite  linear  combinations  of  some  shifts 
of  (j),  they  have  the  same  regularity.  From  their  definition  the  edge  scaling 
functions  generate  all  the  polynomials  up  to  degree  N  —  1  which  ensure  an 
order  N  approximation  over  ail  the  interval,  and  the  existence  of  N  van¬ 
ishing  moments  for  the  edge  wavelets.  With  these  oscillations  and  enough 
regularity,  these  wavelets  basis  form  an  unconditional  basis  for  the  Holder 
spaces  (^^([0, 1])  [1].  The  fast  wavelet  transform  [10]  which  is  essential  for 
most  numerical  applications  is  preserved  even  near  the  boundary  thanks  to 
the  modified  scaling,  (7)  and  detail,  (9)  relations. 


Our  aim  is  to  adapt  this  construction  to  obtain  wavelet  families  generating 
functional  spaces  with  homogeneous  boundary  conditions.  More  precisely  we 
will  consider  the  following  constraints 

y(CLO)(o)  =  0, 

where  is  the  i-th  derivative  of  /. 

As  will  be  shown,  most  of  the  above  construction  will  be  preserved  as 
well  as  numerical  efficiency. 
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II  MULTI  RESOLUTION  ANALYSIS  WITH 
HOMOGENEOUS  BOUNDARY  CONDI¬ 
TIONS 


This  section  is  devoted  to  the  construction  and  the  properties  of  compactly 
supported  wavelet  satisfying  homogeneous  conditions  of  type  = 

yr(CLl)(i)  =  0. 

The  starting  point  has  been  described  in  the  previous  section  and,  keeping 
the  same  notations,  we  now  assume  that  the  compactly  supported  wavelets  on 
the  line  satisfy  0  <  CLO^CLl  <  N  —  I  and  r  >  max{CL0,  CLl).  Therefore, 
the  spaces  Vj([0, 1])  defined  in  (6)  are  included  in  (^^([O,!])  with  r  >  s  > 
max{CL0.,CLl). 

II,1  Construction 

As  in  the  previous  section,  we  only  focus  on  the  left  edge  a;  =  0. 

According  to  (6),every  function  fj  €  1^([0, 1])  is  written 

N-l  2^-N-l  2^-1 

fj(^)  ~  ^  'i'  ^  (12) 

k=0  k=N  k=23-N 

Moreover,  only  the  left  edge  functions  $£,0  are  non  zero  around  x  =  0  and 
then 

(13) 

A:=0 

Therefore  one  way  to  impose  =  0  is  to  enforce  that  all  the  left  edge 

scaling  functions  satisfy  this  condition.  Following  P.  Auscher  (in  [2]),  this 
constraint  can  be  related  to  a  polynomial  behavior. 

Indeed,  from  the  last  section  we  learned  that  is  a  polynomial  of  de¬ 
gree  A^  —  1  on  the  interval  [0,^],  say  for  example  Pj,k{x)  =  -f  a},*®  + 
...  -b  The  CTO-nth  derivative  of  at  0  is  then  equal  to  afk'^- 

Therefore  \o)  =  0  <=>  =  0.  The  construction  of  scaling 

functions  satisfying  =  0  is  then  equivalent  to  the  construction  of 

edge  functions  such  that  their  restrictions  to  [0,2~^]  as  no  component  on  the 
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monomial 

The  first  step  of  our  algorithm  is  then  to  construct  a  family  of  N  edge 
scaling  functions  ^£,0  =  =  0, iV- 1}  and  one  of  N  edge  wavelets 

^£,0  =  ^  =  O5  —5  N  —  1}  with  the  particularity  that  only  one  scaling 

function  and  one  wavelet  contain  on  their  polynomial  part.  The  second 
step  is  to  remove  the  scaling  function  containing  and  to  modify  the 

corresponding  wavelets.  For  simplicity  we  work  on  the  interval  [0,+oo[  with 
a  zero  dilation  scale  {j  =  0),  omitted  in  the  next  notations.  Moreover,  we 
call  pk{x)  =  al  +  alx  +  ...  +  the  restriction  of  on  [0,1]. 

We  start  with  the  first  N  edge  scaling  functions  (5)  of  section  I.  They  are 
defined  with  the  coefficients  so  that 

k  SN— 2ic— 2 

=  i  -  n)  A:  =  0,  ...,7V -1(14) 

n=0  7i=iV 

(see  [1]  for  the  computation  of  these  coefficients).  The  following  proposition 
tells  us  how  to  modify  ipl  to  eliminate  in  the  polynomials  Pk^k  ^  CLO. 
We  call  <p°k  the  new  functions  and  Pk{x)  =  Z)i=o 

Proposition  ILl  The  family  {(p\,  A:  =  0, ...,  TV  —  1}  defined  by: 

i  Tk  =  Tk  k  =  0,...,CL0  .  . 

1  -  WcLo  k  =  CL0  +  1, ...,  TV  -  1 

with 

^0  1  0  \ 

«  _  O‘CL0,k  +  Z^i=CL0+\ 

n  0 

^CLQ,CL0  ~  ^k,k 

is  such  ak^°  =  0,'ik  ^  CLO. 
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Proof: 

The  existence  of  is  always  ensured  for  k  ^  CLO  since 
[1].  Because  pk  is  a  polynomial  of  degree  k,  there  is  nothing  to  change  for 
k  <  CLO,  and  therefore  as  well  as  pkix)  =  Pki^)  for  k  <  CLO.  Given 

k  >  CLO,  let  us  suppose  that  Wl  <  k,l  ^  CL0,a^^'^  =  0.  From  relation  (14) 
we  obtain  the  following  scaling  relation  for  (fl- 

<fl{x)=  - Afca2.ro.)  <^"(2^) 

+Ak,CLo  <PcLo{‘^x)  +  alcLo+i  <^CLo+i(2a;)  +  ...  +  <fc  <pl{2x) 


+  En=;v^*^  “  ^k<y^CL0,n) 


where 


fc-i 

Ak,CL0  =  (^l,CL0  ~  if^CL0,CL0  ~  ^k,k)  +  XI 

i=CI,0+l 


(16) 


Since  (f>{2x  —  n)| [0,1/2]  =  0  for  n  >  iV,  the  contribution  of  the  third  RHS 
term  of  16  to  ^  is  0.  Moreover,  for  0  <  i  <  A;,  (f°{2x)  =  p-{2x)  are  polynomial 
with  no  component  on  Therefore,  the  contribution  of  to  x'^^^  is 

entirely  due  to  Ak,cLo  and  Ak,cLo  =  0  is  the  condition  we  are  looking  for, 
that  completes  the  proof.  ■ 


For  k>  CLO  the  supports  of  (fl  are  no  longer  staggered  but,  in  compen¬ 
sation  VA:,  0  <  k  <  N  —  l,<,3°|[o,i]  is  still  a  polynomial  of  degree  k.  Therefore 
the  functions  <^2,0  <  k  <  N  —  1  are  independent.  Moreover,  they  are  or¬ 
thonormal  to  the  {4>{x  —  n),  n  >  N}  since  they  are  linear  combinations  of 
the  A;  =  0,...,A^- 1}. 

Following  the  previous  section  we  now  orthogonalize  the  family  {(p\-,0  < 
A:  <  —  1),  keeping  the  “monomial  independence.”  The  only  thing  to  do  is: 

to  exchange  the  place  of  <PcLO  before  starting  the  algorithm  from 

index  0  up  to  N-1.  The  result  is  an  orthonormal  family  of  N  edge  scaling 
functions  k  =  0, ...,  —  1}  with  the  particularity  that  only  Pjv_i ,  the 

restriction  of  on  [0,1],  contains  x^^°.  They  satisfy  a  modified  scaling 
equation: 

N-l  3N-2 

<pI’\x)  =  •£  2  lllJ{2x  -n)  k  =  0 N  -  1(17) 

n=0  n=N 
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with  =  0  ioT  k  =  0,  ...,N  —  2. 

The  construction  of  the  N  edge  scaling  functions  for  the  right  edge  comes 
from  the  same  algorithm  for  the  half  line  ]-oo,0].  The  k  =  2^—N-\- 

1, 2-^  —  1}  are  independent  of  and  only  contains  this  monomial 

on  [0,1]- 

After  a  dilatation  of  2^  for  the  0  and  1  edge  functions  and  adding  the 
2^  —  2N  interior  scaling  functions  4>j^ki  one  gets  therefore  a  new  orthonormal 
basis  of  Vj([0, 1]),  the  space  defined  by  (6).  In  this  family,  only  <p^j]cLQ  (resp. 

<p]lcLi)  contributes  to  (resp.  x^^^)  on  [0, 2“-']  (resp.  [1  -  2“'^,1]). 

To  perform  our  first  step  of  construction  we  now  continue  by  isolating  a 
single  wavelet  containing  x^^  on  [0,1/2]. 

As  in  the  previous  section  N  wavelets  at  each  boundary  should  be  added 
to  the  interior  family  Focusing  again  on  the  left  edge,  we  construct  a 
first  family  following  (8)  as 


■^(2®)  -  S  k  =  0, ..., N  -  (118) 

n=0 

Again,  each  function  is  polynomial  on  the  interval  [0,1/2]. 

However,  since  for  all  A:,  depends  on  all  the  0°  contains  the 

monomial  x'^^  and  are  therefore  not  suitable  for  our  first  step  (we  remind 
that  we  want  to  construct  a  family  of  edge  wavelets  such  that  only  one 
contains  x^^^  on  [0,1/2]).  Still,  from  (17)  and  (18)  we  deduce  a  modified 
detail  equation  for  these  functions  that  writes 

N-l  3JV-2 

=  E  +  E  A°n«2^  -  ")  i  =  0, ....  Af  - 1.  (19) 

n=0  n^N 

The  following  proposition  tell  us  how  to  transform  the  functions  to  reach 
our  first  step. 

Proposition  IL2  The  family  =  {^fc.  A:  =  0, ...,  N  —  1}  given  by: 


I  =  lA? - WV-S,-,  A  =  0,...,W-2  A%-i 

\  '^N-1  —  i’N-l  1^-1, N-l 

is  such  that  only  the  restriction  ofip%_i  to  [0,1/2]  contains  x^^. 


(20) 
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Proof:  In  (19)  the  monomial  is  present  only  in  Writing  the 

details  equation  for  and  canceling  the  coefficient  of  gives  the  desired 
result.  ■ 

As  previously,  we  apply  a  Gram-Schmidt  orthonormalization  that  pre¬ 
serves  the  above  property.  Indeed,  starting  from  ipQ  we  get  an  orthonormal 
wavelets  family  ^ =  0,  —  for  which  only  'tj’E-i  contains 

on  [0,1/2].  These  wavelets  are  defined  using  the  details  equation 

N-l  3W-2 

=  E  +  E  G!X2x-n)  Jr  =  0....,JV-l(21) 

n=0  n=N 

It  only  remains  to  make  this  construction  again  for  the  right  edge  with 
monomial  and  to  expand  all  the  boundary  wavelets  of  a  factor  2^.  To¬ 
gether  with  interior  wavelets  family  ^/,  they  form  an  orthonormal  basis  of 

We  have  now  reached  our  first  step  since  we  have  constructed  a  basis 
of  scaling  functions  for  I^  ([0, 1])  and  a  basis  of  wavelets  for  Wj([0, 1])  such 
that  in  each  family,  only  one  function  has  a  component  on  on  [0,1/2] 
and  only  one  function  has  a  component  on  x^^^  on  [1-1/2]. 


As  announced,  we  now  perform  the  second  step  of  our  construction  by 
removing  the  function  the  left  edge  and  the  corresponding  ones, 

^or  the  right  edge. 

The  last  technical  point  is  the  modification  one  should  make  to  the  wavelet 
space.  We  have  the  following  proposition: 
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with  a  and  b  solutions  of 

AT-l  +  =  0 

(24) 

a^  +  b^=  1 

and  a’,  b’  solutions  of  the  same  set  of  equations  with  coefficients  for  the  right 
edge. 

Then  the  new  family is  an  orthonormal  basis  ofWj{[0,l]), 
the  orthogonal  complement  ofVj{[0, 1])  in  1]).  Moreover  every  scaling 

function  ofVj{[0, 1])  and  every  wavelet  ofWj{[Q,  Ij)  satisfies  the  homogeneous 
boundary  conditions  =  0. 


Proof:  We  prove  only  the  result  for  the  left  edge. 

Let  us  first  recall  that  and  are  respectively  two  basis  func¬ 
tions  of  Vj([0, 1])  and  W)([0,1]),  and  that  Vj([0, 1])  -L  W)([0, 1]).  0°  is 

then  orthogonal  to  all  the  other  basis  functions  of  Vj {[0,1])  and  VLj([0, 1]). 
Moreover,!  |0°|  1^2  =  1  if  and  only  if  a^  +  b'^  =  1.  The  same  argument  holds  for 
0]  and  therefore,  with  the  new  definition  of  ^£,0  and  ^£,1,  ^£,0  U  U 
is  a  family  of  2^  orthonormal  functions. 

Using  the  scaling  (17)  and  detail  (21)  equations  we  get 


N-l 


3N-2 


0°=  E  (»<-!,» +  A,  +  E  + 


n=N 


Taking  into  account  (24),  we  get  that  that  0°  is  independent  of  and 

consequently  belongs  to  V}+i([0, 1]).  Since  the  orthonormal  collection  ^e,o  U 
U  ^E,i  generates  a  closed  subspace  of  t^+i([0, 1]),  orthogonal  to  t^([0, 1]) 
and  of  dimension  2-'  =  dimVj+i{[0,l])  -  dimVj {[0,1]),  it  is  by  definition 
Wj([0, 1])  the  orthonormal  complement  of  t^([0, 1])  in  t^+i([0, 1]). 


Remarks: 

All  these  edge  functions  have  the  same  regularity  as  the  initial  scaling 
function  (f>. 
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Moreover  the  2N  —  2  edges  wavelets  constructed  before  we  remove 
and  conserved  their  N  vanishing  moments.  The  modified 

wavelets  0°  and  0]  belong  to  1])  and  are  therefore  orthonormal  to 

all  the  polynomials  included  in  hj([0, 1]).  But  they  have  no  reason  to  verify 
(0°,a;^^°)  =  (0j,x^^^)  =  0,  since  these  monomials  have  been  excluded  from 

the  edges  of  Vj([0, 1]).  Hence  only  one  vanishing  moment  for  one  wavelet  at 
each  boundary  has  been  lost. 

As  in  the  initial  construction,  the  modified  scaling  and  detail  relations 
insure  that  fast  algorithms  related  to  the  different  basis  projections  are  avail¬ 
able. 

At  this  point  however,  we  don’t  know  exactly  what  kind  of  space  the 
multi  resolution  family  t^([0, 1])  approximate.  This  is  the  purpose  of  the 
next  subsection. 

II.2  Approximation  results 

We  now  check  the  intuitive  result  that  the  wavelet  basis  derived  from  the 
last  construction  is  an  orthonormal  basis  for  suitable  homogeneous  spaces 
on  [0,  l].We  give  a  complete  proof  for  the  Dirichlet  homogeneous  boundary 
conditions  /(O)  =  /(I)  =  0,  i.e,  CLO  =  CLl  —  0,  corresponding  to  !])• 

Using  our  construction  for  CLO  =  CLl  =  0,  we  first  obtain  a  subspace 
V^  ([0, 1])  defined  by  the  orthonormal  basis  ^e,o  U  U  ^e,1i  with  the  partic¬ 
ularity  that  only  two  scaling  functions  are  non  zero  at  the  edges.  It  is  known 
that  under  some  specific  conditions  (see  the  previous  section)  V)([0, 1])  be¬ 
longs  to  a  multi  resolution  analysis  of  the  Sobolev  space  i?^([0, 1]). 

Let  us  simplify  the  notations  and  write  =  (p°  f.  and  Then 

we  have  the  following  result: 

Proposition  IL4  Let  lj([0,l])  be  the  subspace  spanned  by  the  orthonormal 
basis 


Assume  these  scaling  functions  have  enough  regularity  to  involve 

LrUiM)  =  ^‘([0.i]) 

i>io 
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and  that  only  ^ 

U  ^i([0>  1])  =  U  (^([0^  1])  -  span{<filN.i]  <^},2i-jv})  =  ^o([0, 1]) 

i>io  i>io 


Proof:  Take  a  function  /  in  iJo([0, 1]),  and  call  Tlj{f)  and  nj(/)  its  or¬ 
thonormal  projection  onto  Vj([0, 1])  and  1^([0, 1]).  We  have  to  establish  the 
relation 

.lim  ||/-nj(/)||H,  =0  (25) 

where  the  ff^-norm  is  taken  as  ||/|l|^i  =  ||/|||,2  +  llffllia-  Following  the 
density  of  Vj([0, 1])  in  1])  this  is  equivalent  to 

lim  Pi)/)  -  nj(/)||H.  =  0 

J-^+OO 

Now  using  the  orthonormal  basis  of  1^([0, 1])  and  t^([0, 1])  we  have 

W^Af)  ~  nj(/)||Hi  =  \{f:'Aj,N-l)\  Ibj.AT-l lli?i  +  \{fi‘p],23-N)\  \\^],2i-N\\H^ 

since  the  support  of  the  0  and  1  edges  scaling  functions  do  not  overlap. 
Because  ||v5°iv-il|L2  =  1  and  V^°iv-i  belongs  to  i?n[0>  1])  to  the  regular¬ 
ity  of  the  initial  function  f),  we  have 

PX/)  -  n,(/)||».  <  Ci2'(|{/, ,>?„_>)!  +  K/.v>,\2, _;,)!) 

where  Ci  is  a  constant  independent  of  j.  Therefore  we  have  to  check  that 
lim  2>|{/,vVi>l=  2i{/.v;,2r-w)l  =  0 

j->-+oo  ^->4-00 

To  see  this,  we  use  the  inequality 

2’l(/.v"«-.)l<ft||/|i„.  (26) 

with  C2  independent  of  j,  which  will  be  justified  at  the  end  of  the  proof.  Let 
us  take  now  a  sequence  of  functions  {fn)neN  convergent  to  /  with  respect  to 


the  if^-norm,  for  example  /„(x)  =  /(x)x[i,i_i].  Applying  the  last  inequality 
to  /  -  /n  we  get 

2^{f-fn,V>lN-l)\<Cz\\f-fn\\m. 

It  only  remains  to  note  that  there  exist  an  integer  J,  dependent  of  n,  for 
which 

Vi  >  J,  2^|(/  - 

Indeed,  for  fixed  n  we  take  J  such  that  the  supports  of  fn  <1°  i^^t 

overlap. 

Making  n  tends  to  +oo,  and  consequently  i,  leads  to 

,lim  2^|{/,¥j;_jv-i)I  =  0. 

j^+oo 

Obviously  the  same  arguments  holds  for  the  scalar  product  2'^i(/,  v?°2>-iv)l 
and  the  proposition  is  proved.  ■ 

We  still  have  to  establish  the  inequality  (26).  An  integration  by  parts 
implies  that 

where  is  a  primitive  of  Since  belongs  to  L^{[0, 1])  we  deduce 

<  C42-A 

This  assumption  and  the  definition  of  the  il'^-norm  leads  to  the  desired  result. 
Thanks  to  this  proposition  and  the  definition  of  IT,([0, 1])  (see  prop  II. .3)  we 
deduce  a  decomposition  of  1])  wavelet  basis, 

^i(lO,ll)  =  V'»([O,l))0H',([O,ll).  (27) 

i>3o 

Remarks: 

The  proof  for  Neumann  homogeneous  conditions  is  similar  and  involves 
the  if^-norm.  More  regularity  is  therefore  needed  for  the  basis  functions  and 
a  double  integration  by  parts  to  the  inequality  corresponding  to  (26).  In  that 
case,  the  approximated  space  is  the  strict  subspace  of  //'^([0, 1])  defined  as 
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Some  mixed  homogeneous  boundary  conditions,  for  example  /(O)  =  /(I)  = 
/(i)(0)  =  =  0,  could  also  be  addressed  with  a  similar  construction.  In 

this  case,  two  scaling  functions  at  each  edge  are  removed  from  1^  ([0, 1])  and 
are  employed  to  modify  the  wavelets  of  lfi([0, 1]).  This  construction  leads 
to  a  characterization  of  the  functional  space 

=  =  /(!)=“  /“’(O)  = /‘'>(9  =  0}. 

Since  the  left  and  right  basis  functions  do  not  interact  at  scale  j,  different 
conditions  could  also  be  taken  at  0  and  1. 

The  following  section  is  related  to  the  numerical  estimates  related  to  our 
construction  and  to  various  topics  connected  to  its  application  for  partial 
differential  equation  problems. 

Ill  NUMERICAL  ESTIMATES 

This  section  is  devoted  to  the  numerical  estimates  related  to  our  construction 
for  two  cases  of  homogeneous  boundary  conditions,  i.e,  the  Dirichlet  condi¬ 
tions  and  the  Neumann  conditions.  All  the  following  computations  have 
been  carried  out  beginning  with  the  initial  compactly  supported  function  (j) 
closest  to  linear  phase  constructed  by  I.  Daubechies  [6]  with  iV  =  4.  Since 
no  explicit  analytic  expressions  exist,  this  function  is  defined  through  the 
filter  coefficients  used  in  the  scaling  equation  (2).  These  coefficients  are 
provided  in  [6]: 

h.3  =  -.07576571478950,  h.2  =  -.2963552764600,  /i-i  =  .4976186676328 

/io  =  .8037387518051,  =  .29785779560531,  /i2  = -.099219.5435766 

hs  =  -.01260396726203,  =  .03222310060405. 

The  corresponding  interior  wavelet  V’  is  defined  using  the  coefficients  gn  of 
the  details  equation  (4)  with  Qn  =  (— l)"^2iv+i-n. 


III.l  Dirichlet  Boundary  conditions 

The  application  of  the  last  section  algorithm  with  CTO  =  CTl  =  0  (Dirich¬ 
let  condition)  leads  to  a  multi  resolution  analysis  of  Hq{[0,  1]).  Three  scaling 
functions  and  four  wavelets  have  to  be  added  at  each  boundary  (see  section 
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II).  These  scaling  functions  are  solutions  of  a  modified  scaling  equation  (17) 
and  are  therefore  characterized  by  the  coefficients  and  The  cor¬ 

responding  numerical  estimates  (computed  on  a  16  decimal  digits  computer 
with  an  error  smaller  than  10“^^)  are  listed  in  Table  1.  The  coefficients  „ 
and  Gl„  which  occur  in  the  modified  detail  equation  (21)  are  listed  in  Ta¬ 
ble  2  and  define  completely  the  edges  wavelets.  All  the  following  figures  are 
obtained  using  the  cascade  algorithm  [6].  The  three  x  =  0  edge  scaling  func¬ 
tions,  as  well  as  the  three  x  =  1  edge  scaling  functions,  are  represented  on 
Figure  2  at  scale  j  =  0.  The  corresponding  wavelets  are  plotted  on  Figure  3. 
Notice  that  due  to  the  lack  of  symmetry  of  the  initial  scaling  functions  and 
wavelets,  the  x  =  1  edge  functions  can  not  be  deduced  from  the  x  =  0  edge 
functions  using  a  simple  transformation. 


111.2  Neumann  Boundary  conditions 

The  same  numerical  estimates  corresponding  to  the  Neumann  conditions, 
i.e  CL0=CL1=1,  are  listed  in  Tables  3  and  4.  The  Figures  4  and  5  repre¬ 
sent  respectively  the  scaling  functions  and  wavelets  of  this  multi  resolution 
analysis. 

Remarks: 

Some  zero  coefficients  are  provided  in  Tables  2,  3  and  4.  They  are  ex¬ 
pected  as  follows:  for  instance  in  Table  3  J?o,i  —  ^0,2  ~  since  the  scaling 
functions  (^0,0  Neumann  conditions,  is  by  definition  constant  on  the 

interval  [0,1],  it  does  not  depend  on  1  and  2  which  are  respectively  poly¬ 
nomials  of  order  1  and  2  on  [0,1/2];  this  leads  77°, 1  —  -^^0,2  =  0  in  Table  3. 
Others  zeros  are  expected  using  the  same  arguments. 

111.3  Quadrature  formula 

In  order  to  use  these  wavelets  basis  for  numerical  purposes  one  question 
needs  still  to  be  answered.  Given  a  function  /,  how  can  we  define  a  projection 
Vj([0, 1]),  i.e,  how  can  we  estimate  a  set  of  coefficients  Cj^k  occuring  in  relation 
(12)  and  corresponding  to  /?  The  solution  proposed  here  aims  to  compute 
an  approximation  of  the  orthogonal  projection  of  /  on  V}([0, 1])  defining 
quadrature  formula  to  estimate  the  coefficients  cj^k  =  /  /<7’o,fc-  We  define 
below  a  quadrature  formula  of  order  TV  —  1  in  the  same  philosophy  as  G. 
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Beylkin  et  al.  ([8])  or  W.  Sweldens  et  al.  ([13]).  We  are  therefore  looking  for 
weight  coefficients  oJi^k  such  that 

fVQ,k  ~  H  (28) 

iz=0 

where  the  {a^,  i  =  0, -  1}  are  iV  given  points  taken  in  the  support  of 
(Pq  ).  and  such  that  the  approximation  is  exact  for  the  polynomials  of  degree 
less  than  or  equal  to  iV  —  1. 

It  appears  that  the  weight  coefficients  toi^k  are  the  solution  of  the  following 
linear  system: 

/ =  E  "..‘W'  I  =  0,...,N-1  (29) 

i=0 

Hence  we  need  to  evaluate  the  first  N  moments  of  every  edge  scaling  function. 

Multiplying  the  modified  scaling  equation  (17)  by  leads  to  the  N  —  1 
equations: 


N-2  .  3N-2  . 

I  X^plki^)  =  Kn  I  ^VS,n(2x)  +  Kn  I  -  Tl)  k  =  0,  2 

n=0  ^  n=zN 


Since  the  moments  of  order  /  of  the  interior  function  <f>,  f  x^<^(x),  can  be  esti¬ 
mated  using  the  classical  recurrence  relation  given  in  [8],  (30)  finally  leads  to 
the  following  linear  system  AXi  =  bi  where  the  N  dimensional  vectors  Xi  and 
Xi{k)  =  fx‘(plk(x) 
bi  are  defined  as  and 

and  where  the  entries  of  the  matrix  A  depend  only  on  the  Hk^^.  We  easily 
checked  that  this  matrix  is  always  nonsingular,  (to  see  this,  use  the  fact  that 
EL'o  KJ  <  1  since  =  1. 

We  first  provide  the  numerical  values  of  the  moments  of  order  /,  0  <  /  <  3 
for  N=4: 

Mo  =  l.OOOOOOOOOOOe+00  Mi  =  -1.45.319345240e-02 
M2  =  2.11177120898e-04  M3  =  4.34510522842e-02 
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Then,  the  entries  of  Xi  for  I  =  0,  —  1  for  the  0  and  1  edge  scaling 

functions  corresponding  to  Figure  1  (Dirichlet  boundary  conditions)  are  listed 
in  Table  5  and  6. 

Using  the  values  of  these  moments  and  the  N  given  points  Cj,  we  find  the 
weights  LJi^k  for  every  edge  scaling  function  solving  the  linear  Vandermonde 
system  (29). 

Remark: 

A  quadrature  formula  of  same  order  has  to  be  used  to  estimate  the  in¬ 
terior  scaling  coefficients  Cj^k  =  f  <k<  2^  —  N— 1  to  preserve  a 

constant  order  of  accuracy  all  over  the  interval. 
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IV  CONCLUSION 

Compactly  supported  wavelets  satisfying  homogeneous  boundary  conditions 
on  [0, 1]  have  been  constructed.  All  the  tools  required  for  the  use  of  these 
functions  for  numerical  approximation  of  partial  differential  problems  have 
been  detailed. 

Even  if  all  this  construction  extends  by  tensor  product  arguments  to 
higher  dimensions,  efficient  handling  of  general  open  sets  with  boundary 
conditions  is  still  an  open  problem. 
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LIST  OF  SYMBOLS  USED 


The  I^TgX  code  of  the  mathematical  symbols  used  is  given  to  clarify  their 
identity: 

•  H''1_0([0,1])  and  a  similar  code  for  L'^{M)  ;  1])  ; 

i='([0,lI);C'((0,ll); 

i^([0, 1]) ;  \tilde{V}_jC[0,l])  and  same  for  H(i([0,  Ij) ;  V} ;  ;  V^([0, 1]) 

;  W'rdo,  11) 

y  Vj  :  \overline{\bigciip_{j  \jn  \Z}  V_j}  and  same  for  f)  Vj,  U 

jez  jez  j  tnZ 

^("l(2A:7r)  :  {\hat  \phi}~-((n)}(2k\pi)  ; 

<f)j^k  ■  \phi_-Cj  ,k}  ; 

:  \phi_{E,0}  ; 

'■  \'tilde{\varphi}'‘{0,\perp}_{j  ,k}  ;  The  same  expressions  are  us¬ 
ing  substituting  \varphi  by  \psi  and  \Phi  by  \Psi. 

Q)  :  \Bigl(''{n}_{k}\Bigl)  ;  (,)  :  Mangle  ,  Mangle;  0 

:  \emptyset 

hl^  :  h~0_{k,n}  ;  hl  .^  ;  and  the  same  expressions  with  uppercase 

HandG. 

y(CLO)  .  f-{(CL0)}  ;  ;  pj^k  ■  P-{j  ,k}  ;  4=^  : 

\Longleftrightarrow 

Xk  :  Mambda.k  ;  gk  ■  \niu_k  ;  (^k,n  ■  Xalpha-O.fkjn}. 

I^k,n-  \beta''0_{k,n}  ;  0°  :  \Tiieta_ j "'0  ;  IMIi/i  •  \IAL{H"1}. 

nj(/)  :  \Pi_j(f)  ;  If:  \frac{\partial  f}{\partial  x} 

Xrii-il-  \clii_{[\frac{lKn},l-\frac{l}{n}]}  ;  :  \Xi_j''0. 

:  \  approx  ; 
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FIG.  1.  The  edge  scaling  functions  Q  and  j  for  N-2.  Q  is  constant 
on  [0;1]  and  (/?q  j  is  a  polynomial  of  order  1. 


Table  1:  The  left  and  right  filter  coefficients,  and  Hl.^,  for 
the  construction  with  Dirichlet  homogeneous  boundary  conditions 


/(0)  =  /(l)=0. 

n 

"k=0  0“ 

1 
2 

4 

5 

6 

7 

8 

9 

10 

k=l  0 
1 
2 

4 

5 

6 

7 

8 

9 

10 

k=2  0 

1 
2 

4 

5 

6 

7 

8 

9 

10 


6.3091928199e-01 

9.65141 14329e-02 

3.9878718285e-02 

-3.9847422688e-01 

-4.6171223341e-01 

-4.3450617898e-01 

-1.6823324724e-01 

3.9512280064e-02 

7.4334018686e-03 

-1.9004116026e-02 

-7.1465107925e-01 

1.0487856031e-01 

-1.0275008552e-01 

-2.2668171853e-02 

-3.4790471529e-01 

-5.4438357904e-01 

-2.1013565982e-01 

7.1146281575e-02 

8.8701254411e-03 

-2.2677220474e-02 

-2.2362773506e-01 

5.0082651220e-01 

-1.1707940876e-01 

-7.0092772722e-01 

-4.8072499606e-02 

4.1374409814e-01 

1.2755670687e-01 

-6.4838220288e-02 

-4.9723065158e-03 

1.27121191106-02 


HL 

5.9981737601e-01 

1.0261005527e-01 

2.2276045032e-02 

-3.8101599048e-01 

-3.5419483844e-01 

-3.7772113863e-01 

-4.0626758997e-01 

-2.2139835764e-01 

1.5553422358e-02 

3.9763630210e-02 

-6.7917200335e-01 

7.0564886681e-02 

-6.1435155066e-02 

1.0664567825e-01 

-6.9557725550e-02 

-3.1266913055e-01 

-5.5031934226e-01 

-3.3164862537e-01 

2.0464047232e-02 

5.2318055022e-02 

-3.4534868269e-01 

3.9554860545e-01 

-5.1663829154e-02 

-3.4828534215e-01 

-5.6934808461e-01 

-1.6719040087e-01 

4.0449494649e-01 

2.8819962720e-01 

-1.4195865812e-02 

-3.6292922913e-02 


Note:  For  the  right  edge  {x  =  1)  the  coefficients  {^o.nl  related  to  the  scaling  function 
and  are  listed  from  right  to  left.  The  case  n=3  corresponds  to  the  scaling  function 
we  have  removed  and  will  therefore  not  appear. 
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Table  2:  The  left  and  right  wavelet  filter  coefficients,  „  eind  Gl  .^^ 
for  the  construction  with  Dirichlet  homogeneous  boundary  condi¬ 
tions  /(O)  =  /(I)  =  0. _ 


n 

Gk,n 

o 

II 

0 

2.01372233346-01 

2.43728752326-01 

1 

5.42297071696-01 

4.58019738396-01 

2 

-7.06456260676-01 

-3.43170832876-01 

4 

3.51083134916-01 

7.43129442866-01 

5 

1.53699735686-01 

-1.47841134716-01 

6 

-1.02591674556-01 

-1.88985031606-01 

7 

-7.56154289616-02 

4.53878712146-02 

8 

5.57249198776-02 

3.39978239186-02 

9 

2.61854704226-03 

-1.56132174596-03 

10 

-6.69453538106-03 

-3.99165014066-03 

k=l 

0 

-6.68484689456-14 

8.02324168766-14 

1 

2.82336928826-02 

3.87452427586-01 

2 

-5.73227755986-02 

-6.76803923436-01 

4 

5.93855368746-02 

-3.95528874706-01 

5 

-1.14746071166-01 

3.70580260476-01 

6 

-2.90555571956-01 

2.81730011406-01 

7 

8.00740103146-01 

-1.04861072956-01 

8 

-4.97049177256-01 

-8.71045292746-02 

9 

-2.95004668496-02 

3.44429846186-03 

10 

7.54204205196-02 

8.80563822006-03 

k=2 

0 

6.20865228366-14 

-5.08038910776-13 

1 

-1.48020849156-14 

-4.01012046906-13 

2 

6.75040528376-02 

1.02436453536-02 

4 

-3.61416278826-01 

1.13095992306-01 

5 

7.86423738986-01 

-5.49537100056-01 

6 

-4.93103067786-01 

7.6956164180e-01 

7 

2.14217134606-04 

-2.89098979396-01 

8 

5.65649582436-02 

-9.03391646676-02 

9 

-1.08692214076-03 

1.23468124506-02 

10 

2.77880771666-03 

3.1565662736e-02 

k=3 

0 

-2.55842586506-02 

-1.65857759906-02 

1 

-6.58767836676-01 

-6.84201439396-01 

2 

-6.83535934546-01 

-6.45852649396-01 

4 

-3.03310121006-01 

2.59936858876-02 

5 

-3.79712891576-02 

-2.78556944186-01 

6 

6.73154404506-02 

-1.52521275856-01 

7 

1.09446702316-02 

8.71628007516-02 

8 

-7.60747579766-03 

7.26384445936-02 

9 

-3.87732156366-04 

-2.85849822616-03 

10 

9.91269813846-04 

-7.30799073036-03 

Note:  For  the  right  edge  (a;  =  1)  the  coefficients  {Go_„}  are  related  to  the  wavelet  V’j  2i-i 
and  are  listed  from  right  to  left.  The  case  n=3  corresponds  to  the  scaling  function  we  have 
removed  and  will  therefore  not  appear. 
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Table  3:  The  left  and  right  filter  coefficients,  and  for 

the  construction  with  Neumann  homogeneous  boundary  conditions 

/(')(o)  =  =  o> _ 

n  Kn  Hln 

k=0  0  7.0710678119e-01  7.0710678119e-01 

1  O.OOOOOOOOOOe+00  O.OOOOOOOOOOe+OO 

2  O.OOOOOOOOOOe+00  O.OOOOOOOOOOe+00 

4  3.7875157560e-01  3.3279641360e-01 

5  4.1933439083e-01  3.1763076592e-01 

6  3.9462541967e-01  3.3872841031e-01 

7  1.5279210977e-01  3.6432796801e-01 

8  -3.5885680933e-02  1.9854306805e-01 

9  -6.7511337557e-03  -1.3947818884e-02 

10  1.7259840309e-02  -3.56587701 13e-02 

k=l  0  -6.3999878989e-01  -6.2399400714e-01 

1  2.7256158164e-01  2.4746950442e-01 

2  2.8255565264e-02  2.1217268985e-02 

4  9.30331 15429e-02  3.1423763181e-02 

5  3.9549596240e-01  1.5138263011e-01 

6  5.4799002074e-01  3.4935923414e-01 

7  2.1165351693e-01  5.4440142107e-01 

8  -6.7366756980e-02  3.2169170065e-01 

9  -9.0159055096e-03  -2.0365603789e-02 

10  2.3049919460e-02  -5.2066376096e-02 

k=2  0  6.3813937293e-02  1.4204323343e-01 

1  -6.2433761707e-01  -5.3003493793e-01 

2  -7.3965386969e-03  1.7695538529e-02 

4  -6.2421804274e-01  -4.0650580195e-01 

5  -1.3488453082e-02  -4.9029535249e-01 

6  4.3699790930e-01  -9.0099784088e-02 

7  1.4222943082e-01  4.3727314163e-01 

8  -6.9374369289e-02  3.0298274581e-01 

9  -5.5998777298e-03  -1.5509344332e-02 

10  1.4316557613e-02  -3.9650941036e-02 

Note:  For  the  right  edge  {x  =  1)  the  coefficients  {H^  ^}  are  related  to  the  scaling  function 
(pj  2J-3  listed  from  right  to  left.  The  case  n=3  corresponds  to  the  scaling  function 

we  have  removed  and  will  therefore  not  appear. 
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Table  4:  The  left  and  right  wavelet  filter  coefficients,  „  and 
for  the  construction  with  Neumann  homogeneous  boundary  condi¬ 
tions  /(^)(0)  =  =  0. _ 


n 

Gk,n 

Gl, 

o 

II 

0 

2.15782819146-01 

2.34949782526-01 

1 

4.95721926486-01 

4.83864987816-01 

2 

7.22384706406-01 

5.28590342426-01 

4 

-3.91855127606-01 

-6.42785063976-01 

5 

-1.28625379316-01 

4.94295970156-02 

6 

9.64575375906-02 

1.22967765676-01 

7 

6.4796296855e-02 

-1.96385185486-02 

8 

-4.71380546276-02 

-1.17824779836-02 

9 

-2.25556189156-03 

7.31274961366-04 

10 

5.76653336506-03 

1.8695658406e-03 

k=l 

0 

5.25175206776-13 

5.58284904496-14 

1 

2.42582076306-02 

2.04580479866-01 

2 

6.03761404386-03 

5.95315163116-01 

4 

2.28480954346-01 

5.40902183996-01 

5 

-5.18100045206-01 

-4.10311831256-01 

6 

5.63585882656-01 

-3.41803127996-01 

7 

-5.26609196446-01 

1.25291830736-01 

8 

2.84109012846-01 

1.00780642426-01 

9 

2.02151691706-02 

-4.17808156416-03 

10 

-5.16817773606-02 

-1.06816163336-02 

k=2 

0 

1.96606561976-12 

-2.72023522866-12 

1 

3.04374601036-12 

-9.03825977366-12 

2 

5.54245945646-02 

1.52132296616-02 

4 

-2.74831942546-01 

1.26070947286-01 

5 

6.08534761576-01 

-5.57440406146-01 

6 

-1.08500263706-01 

7.62874246586-01 

7 

-6.04472331326-01 

-2.86715290906-01 

8 

4.12820474086-01 

-8.84194258366-02 

9 

2.15540560916-02 

1.22673693236-02 

10 

-5.51047541936-02 

3.13625597116-02 

k=3 

0 

1.99416913306-01 

1.87759949376-01 

1 

5.38130128406-01 

6.17934713206-01 

2 

-6.8862102186e-01 

-6.04315612526-01 

4 

-4.20662970066-01 

-3.70193506376-02 

5 

-7.41228790216-02 

-3.84040103496-01 

6 

1.15186836046-01 

-2.00321062816-01 

7 

2.18612861096-02 

1.30948172306-01 

8 

-1.57509239886-02 

1.06914096436-01 

9 

-7.63899130016-04 

-4.33656314416-03 

10 

1.95297226706-03 

-1.10867878006-02 

Note:  For  the  right  edge  (a;  =  1)  the  coefRcients  {Go_„}  are  related  to  the  wavelet  2i_i 
and  are  listed  from  right  to  left.  The  case  n=3  corresponds  to  the  scaling  functions  we 
have  removed  and  will  therefore  not  appear. 
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Table  5:  The  first  four  moments  for  the  0  edge  scaling  function 
of  Figure  2  (i.e  satisfying  Dirichlet  homogeneous  boundary  condi¬ 
tions)^^ _ 


O 

oo 

¥’0,1 

<^0,2 

1=0 

-1.80518411015e+00 

1.58610287947e-01 

1.42577315923e-01 

1=1 

-3.45978990889e+00 

-1.34279194417e+00 

-1.66298580543e-13 

1=2 

-8.12722032481e+00 

-5.41163047176e+00 

7.89305924956e-01 

1=3 

-2.13375831107e+01 

-1.76739256352e+01 

4.70776245989e-h00 

Note:  For  this  case  the  monomial  x  could  be  expanded  as  a  linear  combination  of  9?o,o 
and  ipg  This  explains  the  zero  value  of  the  second  moment  of  2- 


Table  6:  The  first  four  moments  for  the  1  edge  scaling  function 
of  Figure  1  (i.e  satisfying  Dirichlet  homogeneous  boundary  condi¬ 
tions)^ _ _ _ 


V’o.i 

¥’o,2 

1=0 

-2.03632572831e+00 

2.11233673136e-01 

2.34792863276e-01 

1=1 

-5.17867977397e+00 

-1.87776044074e+00 

-4.21625644274e-13 

1=2 

-1.57657032623e+01 

-1.04701396775e+01 

1.27839462595e+00 

1=3 

-5.24899484894e+01 

-4.53961860524e+01 

1.15515969500e+01 

Note:  same  remarks  as  Table  5  for  the  value  of  the  second  moment  of  2- 
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b)The  three  scaling  functions  for  the  right  edge 


FIG.  2.  The  six  edge  scaling  functions  for  the  case  N=4,  represented  at  scale 
j=0  (i.e.  on  [0,+oo[  for  the  left  edge  and  ]—oo,0j  for  the  right  edge).  All 
these  scaling  functions  satisfy  Dirichlet  homogeneous  boundary  conditions  at 


FIG,  3.  The  eight  edge  wavelets  for  the  case  N=4,  represented  at  scale  j=0 
(i.e  on  [0,+oo[for  the  left  edge  and]  -oo,0]  for  the  right  edge).  All  these 
wavelets  satisfy  Dirichlet  homogeneous  boundary  conditions  at  0. 
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